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We discuss the entanglement properties of bipartite states 
with Gaussian Wigner functions. Separability and the pos- 
itivity of the partial transpose are characterized in terms of 
the covariance matrix of the state, and it is shown that for 
systems composed of a single oscillator for Alice and an ar- 
bitrary number for Bob, positivity of the partial transpose 
implies separability. However, this implications fails with two 
oscillators on each side, as we show by a five parameter family 
of explicit counterexamples. 
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I. INTRODUCTION 

Many experiments in the young field of quantum infor- 
mation physics are not carried out on finite dimensional 
quantum systems, for which most of the basic theory has 
been developed, but in the quantum optical setting. In 
that setting the basic variables are quadratures of field 
modes, which satisfy canonical commutation relations, 
and hence have no finite dimensional realizations. It 
would seem that the theory therefore becomes burdened 
with all the technical difficulties of infinite dimensional 
spaces, while theoreticians are on the other hand still 
struggling to answer some simple questions about qubit 
systems. However, the states relevant in quantum optics 
are often of a special kind, and for this class the typical 
questions of quantum information theory are luckily of 
the same complexity as for the usual finite dimensional 
systems. 

This simple class of states of "continuous variable sys- 
tems" is the class of Gaussian states, i.e., those states 
whose Wigner function is a Gaussian on phase space. 
Such a state is therefore completely specified by its mean 
and its covariance matrix, where the mean is irrelevant 
for entanglement questions, because it can be shifted to 
zero by a local unitary (phase space translation). It turns 
out that the basic entanglement properties of a Gaussian 
density matrix (as a state on two infinite dimensional 
Hilbert spaces) can be translated very nicely into prop- 
erties of its covariance matrix (see Section 2), so that 
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problems involving Gaussian states are reduced to prob- 
lems of finite dimensional linear algebra rather reminis- 
cent of the problems involving finite dimensional density 
matrices. 

For the latter it is well known [nU2] that the positivity 
of the partial transpose ("ppt") is necessary for sepa- 
rability, but sufRcient only for the smallest non-trivial 
systems, namely systems in dimensions 2 (g) 2 and 2 (g) 3. 
In all higher dimensions we can find "bound entangled 
states", which are not separable, but nevertheless have 
a positive partial transpose, and are hence not distilli- 
blc 0. In the case of continuous variable systems the 
first nontrivial examples of this kind were obtained in 
Q. In the Gaussian setting it was shown by Simon ||] 
that for bipartite systems with one canonical degree of 
freedom on each side (Alice and Bob), i.e., once again 
for the simplest possible systems, the equivalence of ppt 
and separability also holds. For this system it was also 
shown that non-ppt states are indeed distillible g]. In 
the present paper we settle the relationship between sep- 
arability and ppt for all higher dimensions, showing that 
the equivalence holds also for systems of 1 x A'' oscilla- 
tors, but fails for all higher dimensions. We show this by 
giving explicit examples for 2 x 2 oscillators. 

The key idea for constructing bound entangled Gaus- 
sian states is the notion of "minimally ppt" covariance 
matrices. These are defined as the covariance matrices 
of ppt Gaussian states, which are not larger (in matrix 
ordering) than the covariance matrix of any other ppt 
Gaussian state. It is easy to see that a minimally ppt 
covariance matrix belongs to a separable state iff that 
state is a product state. Hence bound entangled Gaus- 
sians arise from all minimally ppt covariance matrices, 
which are not block diagonal. Numerically, minimally 
ppt covariance matrices can be obtained very efficiently 
by successively subtracting rank one operators from a 
given covariance matrix. This algorithm is reminiscent 
of techniques for density matrices in the context of "best 
separable approximation" [[7|. Running this procedure 
for 2 X 2 or larger systems generically gives bound entan- 
gled Gaussian states. 

Our paper is organized as follows: In Section 2 we will 
set up the basic notation, and the translation of sepa- 
rability and ppt conditions into properties of covariance 
matrices (for separability this appears to be new). We 
also describe the minimally ppt covariance matrices. In 
Section 3 we prove the equivalence for the 1 x A^ case, 
and in Section 4 we present a five parameter family of 
2x2 bound entangled states. 



II. GAUSSIAN STATES AND ENTANGLEMENT 

A system of / canonical degrees of freedom is described 
classically in a phase space, which is a 2 /-dimensional 
real vector space X. The canonical structure is given 
by a 2/ X 2/ matrix a, known as the symplectic matrix, 
which is antisymmetric and non-singular. With a suit- 
able choice of coordinates ("canonical coordinates"), it 
can be brought into a standard form: The 2/ variables 
are then grouped into / canonical pairs (e.g., position 
and momentum), for each of which the symplectic ma- 

, n ) J and all other matrix 



trix takes the form a = 



elements vanish. 

The symplectic matrix also governs the canonical com- 
mutation relations for the corresponding quantum sys- 
tem: if Ra, a — 1, . . . , 2/ are canonical operators (for 
canonical coordinates these are naturally grouped into / 
standard position operators and / standard momentum 
operators), the commutation relations read 



i[Ra,Rl3\ — o-Q^l . 



(1) 



These relations may be exponentiated to the Weyl rela- 
tions involving unitaries VF(C) = exp(i^ ■ a ■ R), where 
^ & X, and £, ■ a ■ R = J^ap^aO'apR/s- These Weyl op- 
erators implement the phase space translations. We will 
assume that they act irreducibly on the given Hilbert 
space, i.e., that there are no further degrees of freedom. 
Then by von Neumann's uniqueness Theorem S the Ra 
are unitarily equivalent to the usual position and momen- 
tum operators in the C^ space over position space. 

For a general density operator p we define the mean as 
the vector m^ — tT{pRa), and the covariance matrix 7 

by 



7a/3 + *0-a/3 = 2 tr p{Ra ~ mal){Rf3 - TUfsl] 



(2) 



which is well-defined whenever all of the unbounded pos- 
itive operators R^ have finite expectations in p. Due 
to the canonical commutation relations the antisymmet- 
ric part of the right hand side is indeed the symplec- 
tic matrix, independently of the state p. The state- 
dependent covariance matrix 7 is therefore real and sym- 
metric. Moreover, 7 -I- icr is obviously positive definite. 

A Gaussian state is best defined in terms of its char- 
acteristic function, which for a general state is f 1— *■ 
tT{pW{^)). This should be seen as the quantum Fourier 
transform |0| of p, and is indeed the Fourier transform of 
the Wigner function of p. Hence we call p Gaussian, if 
its characteristic function is of the form 



1 



tT[pWiO] =exp(im^e-TC^7C) 



(3) 



Here the coefficients were chosen such that 7 and m are 
indeed covariance and mean of p, as is readily verified 
by differentiation. The necessary condition j -\- ia > 0, 
which is equivalent to 7 — icr > by complex conjugation. 



is also sufficient for Equation ^ to define a positive op- 
erator p. We note for later use that a Gaussian state is 
pure iff (cr~^7)^ = — 1 [|l^, which is equivalent to 7 -I- iu 
having maximal number of null eigenvectors, i.e., the null 
space Af — {<i>| (7-1-10-)$ = 0} has dimension (dimAr)/2. 
Note that this null space must always be considered as a 
complex linear subspace of C^-* , the complexification of 
X. For such a complex subspace we denote by 'ShM the 
subspace of X consisting of all real parts of vectors in Af. 
Then a Gaussian state is pure iff ^Af — X. 

Let us now consider bipartite systems. The phase 
space is then split into two phase spaces X — Xa ® Xb , 
where A stands for Alice and B for Bob. This is a "sym- 
plectic direct sum", which means that a — aA ® (tb 
is block diagonal with respect to this decomposition. 
In other words, Alice's canonical operators Ra com- 
mute with all of Bob's. The Weyl operators are nat- 
urally identified with tensor products: ^^(^yi © ^5) = 
^(Ca) ® W{^b)- We call this an /^ x /b system, if 
dimXA = 2/^, and dimXs = 2/^. 

It is clear from (g) and (^) that the covariance matrix 
of a product state is block diagonal and, conversely, a 
Gaussian state with block diagonal 7 is a product state. 
Separability is characterized as follows: 

Proposition 1 Let 7 be the covariance matrix of a sep- 
arable state with finite second moments. Then there are 
covariance matrices ja cind jb such that 



7 



> 









7b 



(4) 



Conversely, if this condition is satisfied, the Gaussian 
state with covariance 7 is separable. 

In order to show the first statement suppose the given 
state is decomposed into product states with covariance 
7'' and mean ni'' with convex weight X^. Then nia = 
X]fe ^kiTT-a Eind, similarly, for the second moments we have 

7a/3 -\- 2mamf3 = ^ Afc (7^^ + 2to^to^) . (5) 

k 

Hence the difference between 7 and the block diagonal 
J2k ^kl^ is the matrix 



iQ/3 



fe 



Xkmlm^fi - 22 \k\emamfi 



ke 



(6) 



which is positive definite, because ^ ^a^p^ap = 
Y.ki ^kMsk - Si)'^ > 0, where Sk = J2a Ca"!^. 

In order to show the converse, let a be the Gaus- 
sian product state with covariance 7^1 © 7b, and let 
7' = 7 — 7A ® 7b > 0. Then 7' is the covariance of 
a classical Gaussian probability distribution P, and the 
characteristic function of the given state p is the product 
of the characteristic function of a and the Fourier trans- 
form of P. Hence p is the convolution of a and P in 
the sense of M, which is the average of the phase space 



translates W{S,)aW{£,)* over ^ with weight P. Since all 
these states will be product states, p is separable. D 



There are different ways of characterizing the partial 
transpose. One simple way is to say that with respect 
to some set of canonical coordinates the momenta in Al- 
ice's system are reversed, while her position coordinates 
and all of Bob's canonical variables are left unchanged. 
In addition, the order of factors in the partial transpose 
of RaRp is reversed when both factors belong to Alice. 
When we replace p in (||) by its partial transpose, we 
therefore find the antisymmetric part of the equation un- 
changed, whereas 7a/3 picks up a factor —1 whenever just 
one of the indices corresponds to one of Alice's momenta. 
Let us call the resulting covariance matrix by 7. Clearly, 
if the partial transpose of p is again a density opera- 
tor, we must have 7 -f icr > 0. But this is equivalent to 
7 -|- icr > 0, where in a the corresponding components 
are reversed, so that a — {—<Ja) ® cfb- This form of the 
condition is even valid if we do not insist on canonical 
variables. Combining it with the positivity condition for 
Gaussian states we get the following characterization: 

Proposition 2 Let 7 be the covariance matrix of a state 
with finite second moments, which has positive partial 
transpose. Then 



"f + ia > , where a = 



-(yA 

(TB 



(7) 



Conversely, if this condition is satisfied, the Gaussian 
state with covariance 7 has positive partial transpose. 

When p is separable. Proposition || shows the existence 
of a block diagonal 7' = 7^ ® 7b with 7 > 7'. Since ja 
and 7b are covariance matrices in their own right, we 
have "/A ± i(JA > 0, and similarly for Bob's side. But 
this means that 7 > 7' > —iu and 7 > 7' > —ia, and 
7 has positive partial transpose, as a separable density 
operator should. We have made this explicit, because it 
shows that it may be interesting to see how much "space" 
there is between 7 and —ia and —ia. This leads to the 
central definition of this paper: 

Definition 1 We say that a real symmetric matrix 7 is 
a ppt- covariance, if j + ia > and j + ia > 0, and that 
it is minimally ppt, if it is a ppt- covariance, and any 
ppt- covariance 7' with 7 > 7' must he equal to 7. 

Note that a minimally ppt matrix 7 is separable if and 
only if it is a direct sum, i.e., if the corresponding state 
factorizes. There is a rather effective criterion for de- 
ciding whether a given ppt-covariance is even minimally 
ppt: First of all, if there were any 7' < 7 with 7' ^ 7, 
we can also choose 7 — 7' = A to be a rank one oper- 
ator, i.e., a matrix of the form A^^ ~ S,aS,p. Then we 
have 7 -|- icr > eA for sufficiently small positive e if and 
only if ^ is in the support of the positive operator 7 + ia. 



The same reasoning applies to ct, so that 7 is minimally 
ppt iff there is no real vector ^, which is in the support 
of both 7 -f i(T and 7 + ia. Rephrasing this in terms of 
the orthogonal complements of the supports, we get the 
following characterization which we will use later on: 

Proposition 3 Let j be a ppt-covariance, and let J\f and 
M denote the null spaces of"f-\-ia and"f-\-ia, respectively. 
Then 7 is minimally ppt if and only if SJfeA/" and ^J\f 
together span X . 

This gives an effective procedure to find a minimally ppt 
7' below a given 7: in each step one subtracts the largest 
admissible multiple of a rank one operator with vector ^ 
orthogonal to the span of dhJ^ and dhJ\f, which is then 
in the supports of 7 -I- icr and 7 -I- ia. In every step this 
will cither increase M or A/", so that a minimally ppt co- 
variance matrix is reached after a finite number of steps. 



III. THE 1 X N CASE 

This section is devoted to the proof that, for Gaus- 
sian states of 1 X A^ systems, ppt implies separability. 
It is clear from the previous section that this is equiva- 
lent to saying that every minimally ppt covariance ma- 
trix is block diagonal, i.e., belongs to a product state. So 
throughout this section we assume that 7 is a minimally 
ppt covariance matrix. 

As a first step we get rid of irrelevant pure state factors 
in the following sense: Suppose that the two null spaces 
have a non-trivial intersection, i.e., there is a $ 7^ with 
^ £ Af nAf. Then (tr - a)^ = {ij - 17)$ = 0, so <i> 
has non-zero components only in Bob's part of the sys- 
tem. So let Xc denote the subspace of Xb spanned by 
real and imaginary part of <&. Then the restriction of the 
state to the subsystem C satisfies the pure state condi- 
tion (its covariance matrix 7c -I- iac has a null vector by 
construction). It follows that the density matrix factor- 
izes: pa,b\c,c = Pa,b\c ® Pc, where pc is a pure state. 
(This conclusion can also be obtained purely on the level 
of covariance matrices, by introducing in Xb a basis of 
canonical variables containing a canonical basis of Xc). 
Clearly, the separability of such a state is equivalent to 
the separability of p, and the covariance matrix restricted 
to Xa ® Xb\c is again minimally ppt. Hence we have 
reduced the problem to the analogous one for the smaller 
space Xa ® Xb\c ■ 

We may therefore assume without loss of generality 
that the null spaces Af and J\f have trivial intersection. 
This means that we proceed by contradiction, since we 
want to prove ultimately that the state is a product of 
"irrelevant pure state factors" . 

Now let 7^ $ e TV and 7^ $ e TV. Then because 7 
is hermitian, we have ($,7$) = (7$,$). Using the null 
space conditions and the skew hermiticity of a, we can 
rewrite this as 



($,(cr-CT)$) =0. 



(8) 



Now the vector {a — cr)<i> must be nonzero, since other- 
wise we would have $ G A/" fl A/". This is a condition on 
the Xy^-components $a of $, since a and a differ only 
on that two-dimensional subspace. By the same token 
the Xyi-component ^a of $ must be non-zero. Hence all 
vectors ((7 — 5)$ lie in the one dimensional subspace of 
CXa orthogonal to ^a t^ 0. The proportionality con- 
stant is thus a linear functional on J\f vanishing only for 
<I> = 0, which means that J\f must be one dimensional. 
By symmetry dim A/" = 1. By Proposition y the spaces 
dhJ\f and dhJ\f together span X, and since they are two 
dimensional, it follows that dimX < 4, i.e., we can have 
at most a 1 X 1 system. For such systems our claim has 
been shown by Simon [p[, and is hence proved. 



IV. 2 X 2 BOUND ENTANGLED STATES 

It was already mentioned in the introduction that nu- 
merical examples of minimally ppt covariances, which do 
not split into ja ffl 7b are easily generated by the sub- 
traction method. In contrast, the subtraction method 
for 1 X iV systems always ends up at a block diagonal 
7. This is rather striking, but not really conclusive, be- 
cause the numerical determination of the null space of 
a matrix which may have small eigenvalues may depend 
critically on rounding errors. We have therefore prepared 
the following all integer 2x2 example 7: 



7 = 
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(9) 



The key to getting simple examples is symmetry, which 
in turn simplifies the verification of the basic proper- 
ties. The most important symmetry in the example is the 
multiplication operator S with diagonal matrix elements 
(1, 1,-1,-1, 1, -1, -1, 1). It satisfies Sa + aS ^ 0, and 
57 — jS. Consequently, j + ia and j — ia — 5* (7 -I- ia)S 
are unitarily equivalent, so it suffices to check the posi- 
tivity and to compute the null space oi^ + ia. We note in 
passing that this unitary equivalence is not necessary for 
bound entangled Gaussians, since generically the spectra 
of 7 -|- «cr and 7 -I- icr are different. 

Further unitaries commuting with the covariance ma- 
trix (p) are the multiplication operator C with diagonal 
matrix elements (1, —1, 1, —1, 1, —1, 1, —1), and the skew 
symmetric operator R with i?i3 = R24 = R75 = Rse = 1, 
and zero remaining entries. All these operators have 
square ±1, and commute with each other and the sym- 
plectic forms up to signs. Therefore, if we start with 



a generic vector fli G A^, the application of R, C, 5' and 
products of these operators yields eight vectors fli, which 
forrn^ a basis of C*. Since these vectors lie in either 
A/", A/" or their complex conjugates we know how 7 acts 
on them and the covariance matrix is thus determined by 
7 = AO^^, where A, 51 denote the matrices consisting of 
column vectors A^ — jClk- The above 7 is generated in 
this manner from 



f^l-(-l,i,2,-3^,l,-^,l,0) 



(10) 



Then the condition of Proposition || is satisfied by con- 
struction, and we only have to verify that 7 + zcr > 0, 
which is again simplified by this operator commuting 
with R. Explicitly, we get the eigenvalues 0, 3 — \/3, 3, 3-|- 
VS, each with multiplicity 2. 

Generalizing this example we can construct a five pa- 
rameter family of bound entangled Gaussian states com- 
muting with R, S and C in the same manner as above. 
We start with a generic vector 

Jli = (-a,z6, c, -id, e, -i/, 1,0) , a, 6, ...,/> 0. (11) 

Then 7 being real and symmetric requires d = (bc+ f)/a 
and from the characteristic function of (fife. A; -|- icrfli) 
we obtain that 7 + id > iff a < ce, where equality is 
ruled out since this would be equivalent to det(f7) — 0. 

States obtained from (11) are all of a non block diago- 
nal form similar to (y), and are hence bound entangled. 
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